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Abstract 

Using the nonrelativistic QCD(NRQCD) factorization formalism, we calculate the 
color-singlet cross sections for exclusive production processes e + + e~ — * J/b + rj c and 
e + + e~ — * J/b + Xcj (J = 0,1,2) at the center-of-mass energy i/s=10.6 GeV. The 
cross sections are estimated to be 5.5fb, 6.7fb, l.lfb, and 1.6fb for rj c ,XcO,Xci and Xc 2 , 
respectively. The calculated J/b + r] c production rate is smaller than the recent Belle 
data by about an order of magnitude, which might indicate the failure of perturbative 
QCD calculation to explain the double-charmonium production data. The complete 
0(a 2 ) color-singlet cross section for e + + e~ — > XcO + cc is calculated. In addition, we 
also evaluate the ratio of exclusive to inclusive production cross sections. The ratio 
of production to J/ipcc production could be consistent with the experimental 

data. 

PACS number(s): 12.40.Nn, 13.85.Ni, 14.40.Gx 

Heavy quarkonium production is interesting in understanding both perturbative and non- 
perturbative quantum chromodynamics (QCD). In recent years the charmonium production 
has been studied in various processes, such as in hadron-hadron collision, electron-proton 
collision, fixed target experiments, B meson decays, as well as Z° decays. Among them, the 
study of charmonium production in e + e _ annihilation is particularly interesting in testing 
the quarkonium production mechanisms, the color-singlet model and the color-octet model 
in the nonrelativistic QCD (NRQCD)[Jj approach. This is not only because of the simpler 
parton structure involved in this process, which may be helpful in reducing the theoretical 
uncertainty, but also because of the spectacular experimental prospect opened up by the 
two B factories with BaBar and Belle, which will allow a fine data analysis for charmonium 
production with more than 10 s e + e - annihilation events in the continuum at a/s=10.6 GeV. 

Recently the Belle Collaboration has reported the observation of prompt J/b via double 
cc production from the e + e _ continuum[2]. For these results, not only the large cross section 
(~0.9 pb) of the inclusive J/b production due to the double cc is puzzling[3], but also the 
exclusive production rate of J/^r] c , <r(e + + —> J/b + Veil)) x ^{.Vc —►> 4 charged) = 

(0.033i^ ± 0.009)pb, may not be consistent with two previous calculations0, 5J, which 
gave a cross section of a few pb for J/b^c- In fact, recent perturbative QCD estimates 
of the J/bcc cross section are only about 0.1 ~ 0.2 pb[B, i?, 0 HJ. So the calculations 
of exclusive cross sections for e + e~ annihilation into J/b^/c and other double-charmonium 
states such as b/bXcj(b = 0,1,2) will be useful to clarify the problem. Experimentally, 
aside from e + + e~ — * J/b + Aco0; Belle ms has also studied processes e + + e —> Xci + A” 
and e + + e~ —> Xc 2 + X, so we hope that the double-charmonium production involving 
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Figure 1: Feynman diagrams for e + e —> J/i/j + T] c (xcj)- 

Xcj(J = 0,1,2) will be detectable in the near future. In the following we will calculate 
the cross sections cr(e + + e~ —> J/-0 + r/ c ) and a(e + + e _ —> J/ijj + y c j) in the leading order 
perturbative QCD. To this order (~ a^) the color-singlet channel is dominant since all color- 
octet channels are of high order of v, which is the relative velocity of the charm quark and 
anti-charm quark in the charmonium, and therefore suppressed relative to the color-singlet 
channel (the relative suppression is at least of order u 4 for the cross sections). Furthermore, 
in order to compare exclusive production with inclusive production rates associated with the 
Xcj charmonium states, we will calculate the cross section cr(e + + e~ —> Xco + cc), and hope 
these ratios will be useful for both inclusive and exclusive production analyses at a/s= 10.6 
GeV. 

We now write down the scattering amplitude in the nonrelativistic approximation to 
describe the creation of two color-singlet cc pairs which subsequently hadronize to two char¬ 
monium states in the e + e“ annihilation process in Fig. Q as DUE] 

A(a + b^QQ( 23 * +1 L Jf )(], 3 ) + QQC s+1 Lj)(j, i )) = Jcr f Jcr L E EE E EE 

LijjzS^jz s l s 2 jk L z Sz s 3«S4 il 

X («§i, «§2 | S'ljjS'ijjz) (L^LuJjz') S'ljjS'ipz | JJipz) ji 3& | l) 

x ( 53 ; S 4 I SS Z ) ( LL Z ; SS Z \ JJ Z ) <3Z; 3* | 1) 

/ A(a + b - Qj(f) + Q k ( f) + Q,(f) + Q t { f)) (L = S), 

{ e* a (L z )A a (a + b - Q,-(f) + Q k { f) + Q,(f) + 4(f)) (L = P), 

( 1 ) 

where (3j;3fc | 1) = 5 jk /PW c , (3Z; 3« | 1) = 5u/PW c , (si;s 2 | S^S^ Z ) > ( s 3 ; ^4 I 
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SS Z ) , (L^L^z] S^S^z | and ( LL Z ;SS Z \ JJ Z ) are respectively the color-SU(3), 

spin-SU(2), and angular momentum Clebsch-Gordan coefficients for QQ pairs projecting 
out appropriate bound states. A(a + b —> Qj(^f) + Qfc(^f) + Qi(^) + Qi(%f)) is the scat¬ 
tering amplitude for double QQ production and A a is the derivative of the amplitude with 
respect to the relative momentum between the quark and anti-quark in the bound state. 
The coefficients Cl^ and Cl can be related to the radial wave function of the bound states 
or its derivative with respect to the relative spacing as 

c, = i-|fl,(0)i 2 , c p = T|fl;( 0 )i 2 . ( 2 ) 

We introduce the spin projection operators Pss z (.P,q) as [TTl fT2 j 

Pss z (p,q ) = s 2 \SS z )v(^ - q- S!)u(^ + q; s 2 ). (3) 

Expanding Pss z (P, q) in terms of the relative momentum q, we get the projection operators 
and their derivatives, which will be used in our calculation, as follows 


Pi%(P, 0 ) = 


Pis z {p, 0) = f(S z )(ft+2m c ), 

PooiP, 0) = ^=75 0>+ 2m c ), 

■[7° f{S z )U> + 2rn c ) - (ft - 2 m c ) ^S z )~f a ). 


(4) 

(5) 

( 6 ) 


4a/2?7 

Then one can calculate the cross sections for the on-shell quarks in the factorized form 
of NRQCDjTj. The cross section for e + + e~ —> J/^ + rj c process in Fig. |T| is given by 


<x(a(>i) + b(p 2 ) -»■ J/VKPs) +Vc(pQ) = 


27ra 2 a 2 |P s (0)| 4 y / s — 16m 2 f 1 


81 m 2 s 3 / 2 


/-i 


where 9 is the scattering angle between p[ and p^, \M\ 2 is as follows 

\M\ 2 = 


2 16384m 2 (t 2 + u 2 -32m 4 ) 


The Mandelstam variables are defined as 

s = (pi +P 2 ) 2 , 

t = (p 3 — pi) 2 = Am 2 — ^(1 — \jl — 16ml/s cos 9), 

u = (P 3 ~ P 2) 2 = 4m 2 - |(1 + \J 1 — 16m 2 /s cos 0). 

The cross section for e + + e~ —> J/iJ) + \ c j process is 

n(a(pi) + b(p 2 ) -> J/^{p 3 ) + XcApa)) = 

27ra 2 a^|P s (0)| 2 |Pp(0)| 2 ^/s — 16m 2 ri 


|M| 2 dcos$, (7) 

( 8 ) 

(9) 
( 10 ) 

( 11 ) 


27m 2 s 3 / 2 


/-1 


|Mj|‘d cos#, 


( 12 ) 
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where \Mj \ 2 for y c0 , Xci and Xc 2 are given by 


|M 0 | 2 = 2048(90112m 40 - 74752m®t - 74752m®w + 23360m°f 2 + 43136m®tw + 23360 m\u 2 
-3152m 4 t 3 - 7600m 4 t 2 w - 7600m 4 fw 2 - 3152m 4 u 3 + 162m 2 t 4 + 444m 2 f 3 u 
+564m 2 f 2 w 2 + 444m 2 7u 3 + 162m 2 w 4 — t 4 u — 3 t 3 u 2 — 3 t 2 u 3 — tu A ) /(3s 7 m 2 ) , (13) 

|Mi | 2 = 32768(1792m 3 + 256m®t + 256 m\u - 56m 4 f 2 - 64 m\tu - 56 m\u 2 - 4 m 2 c t 3 

—20 m 2 c t 2 u — 20 m 2 c tu 2 — Am 2 u 3 + f 4 + 2 t 3 u + 2t 2 u 2 + 2tu 3 + u A )/s 7 , (14) 


|M 2 | 2 = 4096(145408m 40 - 1024m®i - 1024m®« - 2368m®t 2 - 6400m®tu - 2368m®« 2 
+16m 4 f 3 — 208m A t 2 u — 208m 4 tw 2 + 16m 4 w 3 + 24m 2 t 4 + 72 m 2 c t 3 u + 96m 2 c t 2 u 2 
+72 m 2 tu 3 + 24m 2 w 4 — t A u — 3 t 3 u 2 — 3 t 2 u 3 — tu 4 ) / (3s 7 m 2 c ). (15) 


In the numerical calculations, we choose y/s = 10.6GeV, m c = 1.5GeV, a s = 0.26, 
|i? s (0)| 2 = 0.810GeV 3 and |/?p(0)| 2 = 0.075GeV 5 jI3|. and assume that in the nonrelativistic 
approximation nrij m = m Vc = m XcJ = 2 m c . The numerical result for e + + e~ — > J/i/j + r) c is 

cr(e + + e _ —> J/-0 + rjc) = 5.5fb. (16) 


While the numerical result for the cross section of e + + e~ —> J/ip + r/ c is more than a 
factor of six smaller than the experimental dataj2] (with uncertainties due to the unknown 
decay branching fractions into > 4-charged particles for the rj c ), the calculated ratio of 
cr(e + + e~ —> J/il) + ?/ c )/cr(e + + e“ —> J/V’ + cc) ~ 0.037 might be consistent with the 
experimental result with the choice of cr(e + + e~ — > J/ip + cc) = 148fb obtained by taking 
our input parameters. 

The cross sections for J/i^XcJ production at i/s = 10.6GeV are given as 


cr(e + + e - 

-»• J/'fp + Xco) = 6.7fb, 

(V) 

cr(e + + - 

-»• J/ty + Xci) = 1-lfb, 

(18) 

cr(e + + e~ - 

-*• J/^ + Xc2) = l-6fb. 

(19) 


In Fig. [21 we show the cross sections as functions of the e + e center-of-mass energy y/s, 
and we can see that the cross sections for e + + e~ — > J/-0 + r/ c , e + + e _ —> J/i/j + Xci and 
e + + e _ —> J/ip + Xc 2 decrease rapidly as -y/s increases. But the one with J/ip + Xc 2 decreases 
more slowly than that with J/'ijj + ij c and J/ip + Xci- 

At = 10.6GeV if we choose cr(e + + e~ —> Xci + cc) = 18.lfb and cr(e + + e - —> 
Xc 2 + cc) = 8.4fb which were obtained in the fragmentation approximation in Ref.[2|, then 
we have the ratios 


a(e + + e —> J/V> + Xd) 
a(e+ + e~ -»■ Xci + cc) 
cr(e + + e~ -» J/-0 + y c2 ) 
<j(e+ + e~ —> + c2 + cc) 


0.061, 

(20) 

= 0.19. 

(21) 
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Figure 2: Cross sections for a(e + + e —> J// + r/ c ) (solid line) and a(e + + e —> J// + 
Xcj) (dashed line for J = 1, dotted line for J = 2) plotted against the e + e _ center-of-mass 
energy a/s with 7,—Js/ So and -/so = 10.6GeV. 


n(e+ + e -» J/jp + Xci) 
cr(e + + e _ —> J// + cc) 

n(e + + e~ -> J/ip + x c2 ) 
cr(e+ + e~ —> J// + cc) 


0.007, 

0 . 011 . 


( 22 ) 

(23) 


As for the x c0 inclusive double-charm production the rate was not given in Ref.[S|, we 
calculate cr(e + + e~ —> Xco + cc) in a complete form to the C?(a^) order in perturbative QCD. 
We give the amplitude of the first diagram in Fig. 01 for e + + e~ —> y c o + cc as 


M = C(L z )(lL z - 1S Z \J = 0, J 2 = 0)VCl - C -^h(p 2 ) 7 ^(p 1 )-hz(p c ) 

V 3 s 




X [7“U&7 a O; + 7«Pf s ,7 a O> j (p i ,), 


(24) 


where e c = |e , T a is the SU(3) color matrix, the matrix O M is relevant to the on shell 
amplitude and O' 7 is its derivative with respect to the relative momentum between the 
quarks that form the bound state. We can also express the contributions of other three 
diagrams in a similar way, and our numerical results are obtained with the full contributions 
of these four diagrams. Some useful information of the calculation is given in the Appendix. 

We finally get the cross section for this process 


cr(e + + e —> Xco + cc) = 49fb. 


(25) 
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Figure 3: Feynman diagrams for e + + e —> Xco + cc process. 


Then one has the ratio 


cr(e + + e -> J/ijj + Xdo) 
cr(e+ + e~ -> Xco + cc) 
n(e + + e~ -» J/-0 + y e0 ) 
cr(e+ + e _ —> J//> + cc) 
In Fig. QJ we show cross sections for <j(e + + e —> 


0.14, 

(26) 

0.045. 

(27) 


Xco + cc)(solid line) and cr(e + + e 


J/ip + x c o) (dotted line) plotted against the e + e center-of-mass energy y/s with z=^/s/so 
and = 10.6GeV. One can see the ratio in Eq. decreases drastically as the center- 
of-mass energy increases. This is consistent with the result in Fig. |3 We hope the ratios 
between Eq. m and Eq. (ElD could be tested in the near future. 

In summary, despite of many uncertainties due to the relativistic corrections, the QCD 
radiative corrections, the possible color-octet channel contributions, and the choice of phys¬ 
ical parameters (e.g. the charm quark mass and the strong coupling constant), both the 
inclusive and exclusive double charm production cross sections calculated in perturbative 
QCD turned out to be seriously underestimated as compared with data. Therefore we in¬ 
tend to conclude, as in |3j, that it seems hard to explain the double charm production data 
observed by Belle based on perturbative QCD (including both color-singlet and color-octet 
channels), and possible nonperturbative QCD effects should be considered at y/s = 10.6GeV. 

While we were about to submit our result, there appeared one paper which also con¬ 
sidered exclusive double-charmonium production jl 5|. Those authors took the QED effects 
into account in addition to the QCD effects that we considered. We find our result for the 
exclusive double-charmonium production is consistent with theirs but we also analyzed some 
inclusive processes which were not discussed in Ref. na. 
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Figure 4: cross sections for a(e + + e —► y c0 + cc)(solid line) and a(e + + e —> J/-0 + 
Xco) (dotted line) plotted against the e + e - center-of-mass energy y/s with z=\Js/s 0 and 
y/sQ = 10.6GeV. 
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Appendix 


In this appendix we give the cross section for the e + + e —> Xco + cc process shown in Fig. |31 


da = 


m 


h 4 (pi + P-2 - Pc~ Pc~ P) 


2s(27t) 5 

It is convenient to rewrite the cross section as 

\M \ 2 


d 3 p c d 3 p 5 d 3 p 
2 E c 2 E s 2 E 


da = 


4 1 , d 3 p c d 3 pc d 3 p 4 


2s(2tt) 5 

|M | 2 

2s(2tt) 5 


5 (Pi +P 2 ~ V ~ p)S (v ~ Pc - Pc) 


2 E r 2 Er 2 E 


dp 


4 4 d 3 p c d 3 p 5 d 3 pd 3 r] 2 

o (Pi + » 2 — V ~ V)o (V — p c — Pc) -am., 

u v / fc ^ 2 ^_ 2 ^ 2 ^ r,-. 


(28) 


(29) 


where m 2 = r/ 2 
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The integral over the phase-space of cc is evaluated in the corresponding center-of-mass 
frame 


d 3 p' c d 3 p' 5 


sV 


Pc 


■pQ 


1 

8 m 2 


A 1//2 (m 


2 

rp 


m 2 c ,ml)dQ!, 


where A (a 2 , b 2 , c 2 ) = a 4 + b 4 + c 4 - 2 a 2 6 2 - 2a 2 c 2 - 26 2 c 2 . 

The remaining integration are performed in the e + e _ center-of-mass frame 


d 3 p d 3 i 7 
2E 2 Efj 


5 4 (p 1 +p 2 


v-p) 


1 

8 s 


A 1//2 (s , m 2 , m 2 )dQ, 


where rn p = 2 m c , is the mass of the bound state. 
Finally we have 


(30) 


(31) 


da 


\M\ 2 C f 


64s 2 (27r) 5 m^m, 


-A 1 // 2 (m 2 m 2 ,m 2 )A^ 2 (s,m 2 m 2 )dfl'dfldm 


v 


(32) 


The limit of m n is 

2m c < m v < yfs — rn p . (33) 

To accomplish the integration we use the Lorentz transformation between the two frames 
as L = R 2 R\, where 


Ri = 


1 + ^0 0 


0 

0 


1 

0 


a 

m n 

0 

0 


00 


1 + ^r 

m v 


r 2 = 


(l 0 0 0 ^ 

0 cos 9 0 sin 9 


0 0 


0 


V 0 — sin 6 0 cos 9 ) 

The momenta in the e + e~ center-of-mass frame are 

Pi = (Vs/2,0,0, Vs/2), 

p 2 = (Vs/2,0,0,-Vs/2), 

p c - R 2 R 1 PC, 
p, - RiRip'-c 


(34) 


(35) 


(36) 

(37) 

(38) 


( 39 ) 



( 40 ) 


p = {^Jp 2 + m 2 , \p\ sin 6, 0, \p\ cos 9), 

where p' c and p' 5 are the momenta of c and c in the O' frame, and they are 

p' c = (E' c , \pc\ sin 9' cos 6', \p c '\ sin 9' cos 6', \p c '\ cos 6'), (41) 

p' 5 = (Eg, — \pc\ sin 0'cos 0', — \p c '\ sin 6 1 'cos#', — \p c '\ cos#'). (42) 

In Fig. 01 the lower (non-fragmentation) diagrams give very small contributions (about 3 
percent), so for simplicity here we only write down the expressions for the contribution of 
the upper diagrams and give 

|M| 2 = 2 ( 47F ) 4q;2q ^ ( aa + 2 ab + 66). (43) 

Zj ( 

We define pp x = p.p u pp 2 = p.p 2 , PPz = P-Pc, PP4 = P-Pc , P 13 = Pi-Pc, Pu = Pi-Pc, P 23 = P 2 -Pc , 
P 24 = p 2 -Pc■ We notify aa = bb and 

aa = [4(800?7i). 0 s + 800m®pi4pp 2 + 800m®p 24 ppi + 1440m*pp3S + I60mlpi 3 p24pp 3 
+160?n°pi 4 p 2 3PP3 + 1440m°pi 4 pp 2 PP3 + 1440m°p 2 4PPiPP3 + 900 m°pp 3 s 
+184:m^pi 3 p 2 4ppl + 18Am^pupo3Ppl + 856m A pupp 2 ppl + 856m A p24PPiPP 3 

4 O O O o o o o 

+216m c pp 3 s + 56m^pi3p 24 pp 3 + 56m c pi 4 p 23 pp 3 + I 68 m c pi 4 pp 2 pp 3 
+168m 2 c p24PPiPpl + 13 m 2 c pp\s + 2pi 3 p 24 pp 3 + 2 p 1 A p 2 3 ppt)}l[ 8 m 2 c s 2 (Mm 12 
+I92ml°pp 3 + 240m®pp 3 + 160m®pp 3 + 60m A pp 3 + 12m^pp 3 + ppf)\, (44) 


ab = [4(400mJ°s + 400m^ ]3 pp 2 + 400m®pi 4 pp 2 + A00m s c p 23 ppi + 400/n®p 24 ppi 
+ 4 : 00 mlppipp 2 + A89rrP c pp 3 s + 400m®pp 34 s + 480m*pp 4 s + 8 Qrnppi 3 p 24 pp 3 
+80mlpi 3 p 2 4PP4 + 280m®pi 3 pp 2 pp 3 + 440m®pi 3 pp 2 pp 4 + 80m®pi 4 p 23 pp 3 

-F80m®pi 4 p 23 pp4 + M0m 6 c Pi4PP2PP3 + 280m®pi 4 pp 2 pp 4 + 2807n 6 c p 23 ppipp 3 
+A40mlp2 3 ppipp4 + 440m®p 24 ppipp 3 + 280m®p 24 ppipp 4 + 240m®ppipp 2 pp 3 
—400m°ppipp 2 pp34 + 2A0m^ppipp2PP4 + 140 m c pp 3 s + 240m°pp3pp34S 
+476 m & c pp 3 pp A s + 240m®pp 3 4pp 4 s + 140m®pp4S + 40m^pi 3 p 24 pp 3 + 10Am A pi 3 p 2A pp 3 pp A 
+4:0m A p 13 p 2 4Ppl + 20m A pi 3 pp 2 ppl + 288m^p 1 3 PP 2 PP 3 PP 4 + l20m A p 13 pp 2 ppl 
+4Qrn A p 1 4p2 3 pp 3 + 104m^p 1 4 p2 3 pp 3 pp4 + 40 m A p u p 23 ppl + 120m A p u pp 2 pp 2 3 
+288m A p u pp2Pp 3 pp4 + 2 Qrn 4 c p u pp 2 ppl + ^ m tP23PPiPP 3 + ^88m A p 23 ppiPp 3 pp4 
+120 jti a P23PPiPpI + l20m A p24PPiPP 3 + 288m A p 2 4PPiPP3PP4 + 2tim A p 2 4PPiPpl 
— 240m A pp 1 pp2pp 3 pp 3 4 + l4:Am A ppipp2Pp 3 pp A — 2A0m A ppipp 2 pp 3 4PP4 + 108 m A pp 2 pp A s 
+1 AArri 4 j)p 3 pp 34 pp A s + 1087n A pp 3 ppls + 28m 2 c p l3 p 2 4PP 2 3 PP4 + 28rn 2 c p {3 p 2 4PP3Pp4 
+1 2rri 2 p 1 3 pp 2 pp 3 PP 4 + 72m 2 c pi 3 pp 2 pp 3 ppl + 28m 2 c pi A p2 3 pp 2 3 pp4 + 28m 2 j) XA p2 3 pp 3 pp\ 
+72m 2 c pupp 2 pplpp4 + 12mlpi4pp 2 pp 3 pp 2 4 + I2m 2 c p 23 ppipp 2 3 pp4 + 72m 2 c p 23 ppipp 3 ppl 


9 




+72m 2 c p24ppipplppi + 1 2rnlp24pp 1 PP'sPpI ~ 1 AArn 2 pp ] PP 2 PP 3 PP 34 PP 4 + 1 3m 2 c pplpp\s 
+2pi3P24PPlppl + ‘2pi4P23PPippl)]/[3mls 2 (G‘iml 2 + 96 ml°pp3 + 96m l r Q pp 4 + A8rn H c ppl 
+144m°pp 3 pp 4 + A8rnj)p 4 + 8m c pp 4 + 72 m c pp 3 pp 4 + 72 m c pp 3 pp 4 + 8 m c pp A 
+12rrv t c pp\pp 4 + 36mlpplpp 2 + I2m*pp3pp\ + Qm 2 pp\pp\ + Gm 2 pplpp\ + PP 3 PP 4 )]. (45) 
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